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Abstract 



Interpretation of x-ray absorption near-edge structure (XANES) experiments 
is often done via analyzing the role of particular atoms in the formation of 
specific peaks in the calculated spectrum. Typically, this is achieved by cal- 
culating the spectrum for a series of trial structures where various atoms are 
moved and/or removed. A more quantitative approach is presented here, 
based on comparing the probabilities that a XANES photoelectron of a given 
energy can be found near particular atoms. Such a photoelectron probability 
density can be consistently defined as a sum over squares of wave functions 
which describe participating photoelectron diffraction processes, weighted by 
their normalized cross sections. A fine structure in the energy dependence 
of these probabilities can be extracted and compared to XANES spectrum. 
As an illustration of this novel technique, we analyze the photoelectron prob- 
ability density at the Ti K pre-edge of TiS2 and at the Ti K edge of rutile 
Ti0 2 . 
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I. INTRODUCTION 



X-ray absorption fine structure (XAFS) is being used for studying both the electronic and 
the real structure of solids. For high photoelectron energies (E > 200 eV), the extended x-ray 
absorption fine structure (EXAFS) can be intuitively described in terms of backscattering 
of the excited photoelectron by neighboring atoms and the EXAFS analysis has become a 
standard tool for real structure investigations. On the other hand, x-ray absorption near edge 
structure (XANES) is lagging behind in its applications for structural studies, despite some 
promising applications in selected systemsJiTu The reason for this rests in a more complex 
physics hidden behind XANES, resulting both in a more difficult theoretical treatment 
(multiple-scattering, selfconsistency in potentials, non muffin-tin effects) and in a lack of a 
proper intuitive insight into the formation of XANES peaks. The need for involving XANES 
in structural analysis stems partly from the fact that EXAFS is predominantly sensitive only 
to atomic distances and not to bond angles, and partly from the high signal-to-noise ratio 
for some interesting classes of systems which severely limits the ability to extract EXAFS 
oscillations from their spectra. 

A means to interpret XANES in intuitively plausible terms would, among others, facili- 
tate application of XANES spectroscopy in investigations of both the real and the electronic 
structure. Various procedures were applied in the past with the aim to connect XANES 
spectral features with real structure. Among them, let us mention inspecting the effect of 
adding or removing certain atoms in the test cluster ,i investigating the dependence of the 
height of the pre-peak on the geometry of the nearest neighborhood,! calculating XANES 
by summing over many scattering paths@ or employing "direct inversion" technique for ob- 
taining atomic positions and scattering potential from experimental XANES. I In this study, 
we would like to tackle the problem of interpreting XANES spectra from yet another side, 
namely, we want to explore the probability density of the photoelectrons. Such a procedure 
can be — from a certain viewpoint — considered as an answer to the naive question "where 
the XANES photoelectron really is". Answering such a fundamental question would have 
implications for real structure as well as for electronic structure studies. 

At the beginning of this paper, a general consideration of the problem will be outlined, 
the main goal being to define the task in exact terms. Then we will present few mathemat- 
ical formula which describe the problem and determine how the proper probability density 
ought to be evaluated. Next follows a discussion of some practical aspects of evaluating the 
photoelectron probability density and the technique is illustrated on the pre-edge structure 
of Ti K edge of TiS2 and on the whole Ti K edge XANES of TiC>2. Some more technical 
details are given in the appendices. 

II. LOCALIZATION OF THE PHOTOELECTRON 

XAFS arises due to the energy-dependence of the core-electron photoeffect: An electron 
absorbs an x-ray photon, is ejected off the atom and starts to travel inside the solid. In a 
stationary picture, the wave function of the excited photoelectron can be viewed as being 
subject to multiple scattering by neighboring atoms. The probability of photoabsorption 
oscillates with energy — one can, especially in the EXAFS regime, interpret this oscillatory 
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behavior intuitively as a consequence of either constructive or destructive interference of 
the photoelectron wave function. One can thus associate a particular spectral peak with 
scattering of a photoelectron of a certain energy. The topic of this section (and of this paper 
as a whole) is to explore the spatial localization of this photoelectron. 

A. Wave function problem 

In the framework of the quantum theory, the popular question "where is the electron" 
cannot be answered. However, it is possible to ask what is the probability density P(r) that 
a quantum-mechanical object can be found at a given place r — it is just the square of the 
modulus of its wave function, 

P(r) = Mr)\ 2 . (1) 

So the problem is reduced to the task of finding the wave function ip p h e (r) which would 
appropriately describe the photoelectron participating in the formation of XANES at a 
given energy. 

The correct form of this wave function ought to emerge from the way of calculating the 
x-ray absorption spectrum. The x-ray absorption spectrum (XAS) intensity is proportional 
to the probability w that a photon is absorbed by a core electron. That can be expressed, 
within the first-order perturbation theory, as 

w = 2tt JdvK^lH^l 2 8(E c + uj-E u ) , (2) 

where Hj is the interaction hamiltonian perturbing the initial electron state |0 C ) and the 
sum/integration over v spans any complete set of electron wave functions \ip u ) (Rydberg 
atomic units are used throughout this paper, taking m = 1/2, h = 1, e 2 = 2, c = 2/ at, 
a = 1/137.036). Applying standard procedures, Eq. @ can be transformed into expres- 
sions which involve either sums over Bloch states,^ or molecular cluster basis functions,^ or 
photoelectron diffraction states0 or which, in the case of Green function formalism, dispose 
of the final state wave functions altogether.!!! The important factor is that the choice of the 
set of wave functions \ip u ) does not affect the outcome of Eq. (|2]). This is favorable on the 
one hand, as one does not have to care about the particular form of wave functions \ip u ) 
when calculating XAS intensity. On the other hand, it means that the proper photoelectron 
wave function cannot just be borrowed from Eq. (fj) or from any of its clones. 

A close look at Eq. (^|) reveals that the situation is even worse at the first sight. Namely, it 
follows from Eq. (|2]) that the total absorption rate is actually resultant from many incoherent 
processes. Hence, there is simply no single wave function like iV'phe) which could have been 
inserted into Eq. ([J). Nevertheless, one still can ask what is the probability that any electron 
ejected as a results of an absorption of a photon with energy u can be found at r. One only 
has to reformulate the problem slightly: Instead of searching for |^ p he(^)| 2 , the quantity of 
interest should rather be a weighted sum of probability densities of those wave functions 
which describe states participating in the absorption process,li2l 

F(r) = E^/I^/WI 2 • (3) 
/ 
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The weights Wf with which the participating wave functions contribute to P(r) are the 
probabilities that a core electron |0 C ) is ejected into the state \ipf), 



w f ~ KfclHMf)? . (4) 

They ought to be normalized so that their sum yields the total XAS probability w. Note 
that all this is just another way of saying that the ejected photoelectron is described not by 
a single quantum state but rather by a density matrix, 

v = 5>/|^/)(iM • (5) 

Unfortunately, there appears to be a complication resulting from the use of Eqs. 
Unlike in the case of Eq. (g), the outcome of Eqs. now, namely, does depend on the 

choice of the set of wave functions [complete on the energy surface determined by the 
^-function in Eq. (fj)]. There is no way to fix this choice by relying solely on Eq. (0). It is a 
matter of physical intuition, not formal mathematical procedures, to establish which set of 
states 1-0/) is relevant to the physical process in question. A guide for this choice could be 
the conservation of electron number: As one has one electron in the initial state (it is the 
core electron \<p c )), one has to end up with one electron in the final state as well. 

The choice of states \ipf) is quite straightforward in the case of transitions to bound 
states (say, of a molecule). The final states \ipf) are just asymptotically decaying bound 
states \ipb)i normalizable so that / d 3 r|-?/^(r)| 2 = 1. The situation is more complicated for 
transitions to the continuous part of the spectrum, as there is no a priori preference for 
normalization and/or boundary conditions which such a wave function ought to observe. 

In order to find the correct wave functions, let us contemplate a finite cluster of atoms (a 
situation tacitly assumed in most applications of the real-space formalism)0Hll and follow 
the fate of the initially core electron. As a result of absorbing a photon, this electron is 
torn off an atom and, having its energy above the continuum threshold, must finally leave 
the cluster and turn into a plane wave with a well-defined momentum direction k. This 
means that the elementary incoherent processes involved in the sum (H) must be photoelec- 
tron diffraction events. Indeed, x-ray absorption is conceptually nothing else but angularly 
integrated photoelectron diffraction (PED). The final states \ipf) are, therefore, time- 
inversed scattering states, [t/4 ), which are in turn solutions of the Lippman-Schwinger 
equationBSB 

V4~V) = e ifcr + J^r'G[-\ry)V{r')^ k \r') , (6) 

where Gq (r,r') is the advanced free electron Green function and V(r') describes the po- 
tential of the cluster. The states are normalizable to the delta function 5(k) [apart 
from the constant factor l/(27r) 3 / 2 ], which guarantees that they describe exactly one electron 
at a timeS 

The probability density P{r) obtained via Eqs. cannot be normalized to one. 

However, due to Eq. ([]), it can be related to the probability density of a free electron, which 
is described by the wave function exp(ifcr) and holds thus a constant probability density 
everywhere. So the probability density P{r) is actually measured in "units of free-electron 
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probability density". Thus, by pegging the normalization of \ipf) to the normalization of 
the free-electron wave function exp(ifcr), one keeps a universal definition of P(v). 

Intuitively, the wave function iV'jt ) can be viewed as that wave function, from which 
a plane wave evolves within a sufficiently long time. It represents the state into which 
the core electron "jumps" as a result of the electromagnetic perturbation Hj. Thus, the 
quantity P(r) ought to be interpreted as the probability, that the electron ejected from 
a core level can be found at position r, "just after" having absorbed an x-ray photon. 
By evaluating P(r), one provides the most sensible answer to the seemingly naive question 
about the localization of the XANES photoelectron. We bear in mind, at the same time, that 
this Lippman-Schwinger-like description cannot fully substitute for a proper time-dependent 
treatment. 



B. Evaluating photoelectron probability density P{r) 

In this section, we present equations necessary for calculating the photoelectron probabil- 
ity density P(r) . Although some of them can be found in that or other form in various papers 
dealing with x-ray absorption or photoelectron diffraction theory (especially in Refs. |l6| , |l8l) , 



we present them here anyway in order to embed them into the context of this work, to offer 
the reader a complete set of equations which might be helpful for practical calculations and, 
last but not least, to unify various notations and conventions. 

The proper mathematical expression for evaluating the probability density of electrons 
participating in x-ray absorption process can be obtained by inserting the wave function 
|V4 -) ) into Eqs. (g)-®. We get 



da 



0~V AC J 



P(r) = — ld 2 k 



^\r)~ , (7) 



where the PED cross section da/dQk stems from the partial probability Wf of Eq. (|J) and 
the XAS cross section a XA s, 

ensures correct normalization. 

By keeping only the dipole and quadrupole terms in the electromagnetic hamiltonian 
Hf, the PED cross section can be written as 



da 1 

-auk 



\(^ ] \ir\4> c )\ 2 
+ K^KeVX^c)! 2 ! , (9) 



where e is the polarization vector of the incoming radiation and q is its wave vector 
(q = lflaujq). Employing the muffin-tin approximation, the wave functions ipk~\r) can 
be expanded inside the j-th muffin-tin sphere as 

= Y.^\k)nf{kr)Y L {r) , (10) 
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where single-sphere solutions of the radial Schrodinger equation Tl$\kr) are normalized so 
that they smoothly match the free-space solution 



TZ^(kr) = cot S^' ji(kr) — n e (kr) 

outside the muffin-tin sphere. The double-subscript L stands for the pair (£,m). 
Following the formalism of Ref. the coefficients f3$\k) can be expanded as 



00 



(11) 



(12) 



where the amplitudes P^}{L") satisfy 
The scattering matrix W is an inverse matrix to 



(13) 
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LV 



sin^V^ 



Wlu + 47T (1 - <Sy) £ ih<t;\k\RV\) Y Ll (R l >) C L L ' Ll , 

Li 

(14) 



where the Gaunt symbol Cfv stands for 

d 2 nF L (n)F L *(n)F Ll (n) , 
the free-electron propagator J£^, is 

Li 



and R tJ is defined as 



The amplitudes and the scattering matrix W^ L i are "incoming- waves" analogs of 

the amplitudes B$\L") and the scattering matrix [(T + if) -1 ]^, employed in Ref. |16| for 
analyzing the scattering of an electron by a molecule. See Appendix ^ for a more compre- 
hensive comparison. 

Employing the amplitudes (3$\L"), the PED cross-section can be expressed via 



R 13 



R l - R> 



(15) 



(16) 



(17) 



da 
dfta 



££(-if *Y L „{k)D LLc 



L L" 



+ 
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-aw 



L L" 



(18) 



where the dipole and quadrupole matrix elements are 
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and 



D LLc = Jdrr 3 nf\kr)<f) c (r) Jd 2 rY£(r) efY Lc (f) 

Qll c = Jdrr A uf\kr)<P c (r) 

x fd 2 rY£(r) (if) (qf)Y Lc (f) 



(19) 



(20) 



and L c specifies the angular momentum of the core state located at the central atom R°. 
The XAS cross section follows from Eq. (pi) and Eq. ([3~8|) as 



a XAS = 4™fc W EK 0) (^' 

L" I L 



i 2 2 



D LLc 

2 1 



(21) 



resembling in this form analogous expressions presented, e.g., in Refs. [TD| , P~B| (note that we 
ignore the electron spin throughout this paper). Again, in Appendix [A] we will present few 
other equivalent formulations of Eqs. and (|2~1~D . 

Considering Eqs. (0) and (]T0|), one obtains the expression for probability density of 
ejected photoelectron inside the j-th muffin-tin sphere as 



P(r) 



as. 



d 2 k 



da 



J2^\k)TZ^(kr)Y L (f) 



(J), 



(22) 



where the PED and XAS cross sections have to be taken from Eqs. ( |i~8|) and (pl|). 

In electronic structure studies, it might be helpful to have a tool for investigating the 
angular-momentum character of photoelectrons with respect to a site R J . That can be 
achieved by inserting into the sum ([|) only the angular-momentum projected parts of pho- 
toelectron wave functions, 



f 



(23) 



where stands for the relevant projection operator. Analogously to Eq. (|22]) , one gets 
the probability density of finding a XANES photoelectron at position r with an angular 
momentum £ with respect to the site K* as 



PAr) 



1 



OS, 



d 2 k 



da 
dQ k 



Y,^)nf{kr)Y im {f 



(24) 



meaning that the sum Y^ L in Eq. (|22|) was just substituted with Ysm- 



7 



C. Practical aspects 



This study is motivated by an effort to understand x-ray absorption spectra and, in 
particular, to develop a means of connecting spectral and structural features. It is thus 
desirable to explore how P(r) depends on the energy of the photoelectron. A simple way 
to extract the fine structure from this dependence is to subtract from P(r) the probability 
which would correspond to a single isolated atom — just like the EXAFS can be extracted 
from a raw absorption spectrum by subtracting from it the atomic part (sometimes called 
AXAFS or atomic XAFS).0 The single-atom probability density P a to( r ') can be evaluated 
following the procedure outlined in Sec. [11 B| , taking into account only a single scatterer 



(cf. appendix |B|). By investigating the difference P(r) — P a to(f ), one can see more clearly 
the effects of surrounding atoms on the formation of both XANES and the photoelectron 
probability density. 

In practice, one often wants to compare the importance of particular atoms for generating 
XAFS. For that purpose, it is sufficient to compare not P(r) but rather its integrals inside 
suitably chosen spheres. We can define atomic-site-related quantities and APjfo by 



R 0) 







pW = _ / drr 2 jd 2 fP{r) , (25a) 



APjil = W) J Q N dr r 2 J d 2 f [P(r) - P ato (r)] (25b) 

where R$ is a suitably chosen normalization radius and is the volume of the normal- 
ization sphere around the j-th site. Note that in the case of muffin-tin approximation, only 
spherically-averaged values make sense anyway. 

From the XANES analysis point of view, the quantities P^\ AP a ^ contain still quite a lot 
of unnecessary or "ballast" information. This is due to the fact that the probability density 
P(r) defined by (|7|) is dominated by isotropic density of states (DOS) effects, which are not 
specific to the site from which the photoelectron is ejected^ This DOS-like contribution can 
be quantified by defining a DOS-like probability density P DOS (r), which differs from P(r) by 
assuming a fc-independent or "unidirectional" photoelectron diffraction cross section da/dfl 



asH 



Oxas all 



7- 

47T 



^~\r) 2 . (26) 



Thus states with different k contribute to P DOS (r) with identical weights, just as is the case 
of local r-dependent DOS. It can be shown easily that P D os(^) is indeed proportional to the 
density of states n(r, E), 

n(r,E) = --lmG^\r,r;E) = -^P DOS (r) . (27) 
An atomic-sphere related quantity APdos can be defined analogously to Eq. (p5|), 
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1 rR U) r 

*pgL = vu) I N dr r2 J d2f ( p ( r ) - ■ ( 28 ) 

Typically, P does not differ from P DO s by more than 10% (but often much less). Specific 
XAS-related effects may thus be obscured in P(r) by more general DOS effects. The differ- 
ence probability density, P(r) — P D os( r ), informs how the spatial localization of a XANES 
electron differs from the spatial localization of a "generic" electron (with the same energy). 
By investigating the difference probability APno S , one can filter out effects which are not 
specific for XAS. 

Generally, the total probability density P^ informs where the most of the photoelec- 
tron is located, while the difference probability density APdos is especially sensitive to the 
photon polarization e and to the position of the photoabsorbing atom. We postpone a fur- 
ther discussion of these concepts to Sec. A| , where they will be illustrated on a concrete 
example. 



III. APPLICATION OF PHOTOELECTRON PROBABILITY DENSITY 

ANALYSIS 

We explore the potential of photoelectron probability density (PEPD) analysis by ex- 
amining several XANES spectra, which were subdued to investigation of the origin of their 
peaks in the past. In particular, we will concentrate on the pre-peak at Ti K edge XANES 
of TiS2 and on the whole Ti K edge XANES of rutile TiC>2. 

All the calculations presented here were done for a non self-consistent muffin-tin potential 
constructed via Mattheiss prescription (superposition of charge densities of isolated atoms). 
Using non self-consistent potentials is not a serious drawback in this case as our aim is not 
to achieve the best reproduction of experimental spectra but to demonstrate a method how 
the calculated spectra can be analyzed. 

The exchange-correlation potential of Ceperley and Adleri was used for atomic calcula- 
tions of occupied states. In constructing the Mattheiss potential appropriate for unoccupied 
states, an energy-independent Xa potential with the Kohn-Sham value of a = 0.66 was 
used.S Structural data were taken from the crystin database.^ Only dipole transitions 
were taken into account. This limitation is justified in this study because it was demon- 
strated experimentally that quadrupole transitions do not contribute to the TiS2 pre-preak 
intensity significantly^ and in the case of Ti02 spectrum our concern is with the main 
peaks and the extended XANES region, where the quadrupole contribution again is negligi- 
ble. Muffin-tin radii of nonoverlapping spheres were determined so that single-site potentials, 
which were being superimposed, matched at the touching points ( "matching potential con- 
dition"). The muffin-tin zero was set to the average interstitial potential. The influence 
of the core hole left on the central atom by the excited electron was taken into account 
by calculating the central atom charge density with one electron being moved from the Is 
core level to the lowest unoccupied atomic orbital (a relaxed and screen model). A more 
thorough discussion of how the potential is constructed can be found, e.g., in Refs. |5|,|2"o]. 

When evaluating the sphere-averaged quantities P^\ AP^, and APnos according to 
Eqs. fl2"5p and (ESI), one has to choose the normalization radius for each of the spheres. 
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Throughout this section, we always take identical for all atoms of a given compound and 

(i) 

equate it with the smallest of muffin-tin radii. We found that changing affects rather 
the overall magnitude of PEPD than its fine structure. Consequently, the overall picture, 
as presented in Sec. fll AHnTB| , does not depend on the choice of R$. 

Experimental x-ray absorption spectra are broadened with respect to "raw" theoretical 
spectra because of various many-body processes (finite core hole lifetime, extrinsic photo- 
electron losses etc.). It is thus desirable to anticipate this smearing when evaluating the 
PEPD's — otherwise, one would be overburdened with too many details with little or no 
physical significance. Therefore, we modified the scattering potential by adding to it a small 
negative imaginary part. We set its magnitude so that it simulates lorenzian broadening 
equivalent to the half of a natural core hole width (taken from compilation of Al Shamma 
et al)£Q Such a procedure would certainly be insufficient if one tried to describe the spectral 
broadening in a realistic way — apart from too small core hole smearing, energy- dependent 
inelastic energy losses of the photoelectron are unaccounted for altogether. However, as 
our main aim is to analyze the peaks in the calculated spectrum, we prefer to include less 
damping in order to have certain features more pronounced. We checked that increasing the 
imaginary potential twice would not affect the outcome of our analysis considerably. 



A. Pre-peak at Ti K edge of TiS 2 

TiS2 has a layered crystal structure consisting of sulphur-titanium-sulphur slabs. Each of 
these slabs is formed by a twodimensional hexagonal titanium sublattice sandwiched by two 
closely adjacent sulphur hexagonal sublattices. Locally, the titanium atom is octahedrally 
coordinated by six sulphur atoms (cf. figure 3 of Ref. p7[) , the second coordination shell is 
formed by Ti atoms of the same hexagonal sublattice to which the central Ti atom belongs. 

The polarized Ti K edge XANES of TiS2 shows a distinct pre-peak at the xy polariza- 
tion (i.e.. when the polarization vector of the incoming radiation lies within the titanium 
layer) .0@ Relying on XANES calculations for artificial trial structures, Wu et a@ suggested 
that the pre-edge is generated by multiple scattering which involves mainly the central Ti 
atom and the Ti atoms of the second coordination shell. 
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FIGURES 




FIG. 1. Theoretical pre-peak structure of polarized Ti K edge XANES of TiS2 (upper panels), 
together with the photoelectron probability density P^ (lower left panel), atomic difference prob- 
ability density AP a j„ (lower middle panel), and DOS difference probability density AP^os (lower 
right panel) around selected atoms. Solid lines correspond to the e\\xy polarization, dashed lines 
to the e\\z polarization. PEPD curves are identified by the chemical type of the appropriate atom 
and by its distance in A from the center of the 135-atoms cluster. Thin dotted lines mark zeros of 
P^\ AP^il, and AP^qs f° r each subgraph. The absolute scale of the PEPD's is indicated at the 
lower panels, the scale of XANES is arbitrary. Note that curves for Ti atoms were divided by 20, 
as indicated. 

A comprehensive analysis of the Ti if edge XANES of TiS2, including comparison be- 
tween theory and experiment, can be found elsewhere .0110 Therefore, we present in Fig. 
[I] only the theoretical polarized pre-peak structure, calculated for a cluster of 135 atoms, 
together with P^\ AP£, and AP$ S curves for the central atom and for atoms of its four 
nearest coordination shells. Two polarizations correspond to the e vector averaged over the 
full 2tt angle within the xy plane and to the e\\z setup. Titanium probability densities are 
scaled down by a factor of 20 with respect to corresponding sulphur curves, as indicated. 
Individual atoms belonging to the same coordination shell give rise to identical PEPD's in 
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this case (also for further shells than those displayed in Fig. [T]). As we did not perform a 
band-structure calculation, the Fermi level is not fixed. It follows from Fig. [I] that it ought 
to be around ~ 6 eV above the muffin-tin zero, therefore, all states bellow it are actually 
occupied and do not contribute to x-ray absorption spectrum — we show them just for 
completeness. 

The first information one gets from Fig. [3] is that the photoelectron probability density 
clearly oscillates with energy. The fine structure in PEPD differs from the XAFS. One might 
be surprised at first by the fact that this is the case not only for non-central atoms but also 
for the photoabsorbing one (after all, XAS is approximately proportional to unoccupied DOS 
angularly projected on the central atom). One has to recall, however, that PEPD and XAS 
carry in fact different kind of information: While the XANES intensity informs us about 
the probability that the ejected electron goes anywhere, PEPD tells us how this "anywhere" 
looks like. Only a tiny fraction of the photoelectron density matrix (|5|), namely, that part of 
it which overlaps with the core of the central atom and which has the angular-momentum 
character conforming to the dipole and/or quadrupole selection rules, enters indirectly into 
expressions for the XANES intensity ( |A6| )-([A7|). 

One can see immediately from the lower left panel of Fig. [I] that the chemical type governs 
the gross shape of PU> curves (we checked that this is true also for more distant atoms, which 
are not displayed here). The differences between PEPD's around atoms at crystallograph- 
ically equivalent positions (say, Ti or S atoms belonging to different coordination spheres) 
comes first of all from the fact that this quantity is specifically related to the position of 
photoabsorbing atom, i.e., the site from which the photoelectron was ejected. (Apart from 
that, our finite-cluster approach obviously introduces inequality among otherwise equivalent 
sites.) 

Hardly any polarization effect can be noticed in the or AP a ^ curves, despite the 
fact that the e-dependence in the XANES spectrum is quite significant. Only the difference 
probability APd OS , which emphasizes the effect of the particular site from which the pho- 
toelectron is ejected, displays a strong polarization dependence. This is a manifestation of 
the dominance of the DOS-generated unidirectional Poo's contribution to the total PEPD, 
as mentioned in the end of Sec. |11 C\ 

Neither local DOS nor Puos take into account that the photoelectron is ejected from a 
particular site via a dipole transition^! — they are concerned with all electrons of a given 
energy equally. On the other hand, the probability P^ that a Ti if edge XAS photoelectron 
will be found near atom j differs, albeit slightly, from the probability of finding there "any" 
electron of the same energy. It is this small difference which reflects the fact that the 
TiS2 crystal does not look identical when viewed from the Ti site either in the xy-plane 
or in the z-axis directions. So one has to resort to the difference probability AP^os if the 
polarization-related features of PEPD are to be studied. 

On the other hand, simple subtraction of the single-atom probability density P^l does 
not provide a new insight. Due to the smoothness of atomic probability density P^l, the 
P® and APjto curves look very similar, as can be seen in the left and middle panels of Fig. 
[I]. The total probability density P^ and atomic difference probability density AP^o carry 
essentially identical pieces of information. 

As it follows from Fig. |I|, there seems to be no simple correspondence between XANES 
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peaks and peaks in the photoelectron probability density. Although for some atoms and/or 
peaks, one can establish a visual connection between XANES and PEPD, for other features 
such a discernible connection is clearly absent. This lack of simple correspondence between 
XANES and PEPD peaks may be a manifestation of the interference nature of XAFS — it is 
generated not just through accumulating electrons here and there but rather by interference 
between many scattering paths. 

We performed also an angular-momentum analysis of PEPD, according to Eqs. (P3])-(P^). 
We found that the d component dominates at Ti atoms (it contributes by more that 95% 
to either or AP^to or APpos) and that the p component prevails at S atoms (comprises 
60-70% of PEPD). Moreover, practically all of the fine structure in PEPD is formed by the 
dominant ^-component (i.e., d at Ti and p at S atoms). This might again look surprising 
for someone who is accustomed to the conventional slang that the excited photoelectron 
has a p character at the K edge due to the dipole selection rule. The point is that this 
would apply literally only in case that the final photoelectron state were an eigenstate of the 
angular momentum. On the other hand, the wave function |^ T ), which describes the exited 
photoelectron (as argued in Sec. II A ), is not an angular-momentum eigenstate. Rather, it 
is a superposition of states with different angular momenta, and the component with £=2 
dominates at the central Ti site as a whole, while the relatively tiny component with £=1 
determines the XANES intensity. 

As mentioned, Wu et alii suggest that the second-shell Ti atoms play a crucial role in 
generating the e\\xy pre-peak. They arrived at this conclusion by observing that this pre- 
peak disappears if those six Ti atoms are removed from the cluster. A closer look at Fig. [I] 
reveals that the DOS-inclusive probability is, indeed, much higher near Ti atoms than 
near S atoms. However, this effect is clearly DOS-related, having little connection with the 
particular direction of the photoelectron trajectory. By inspecting the DOS-corrected differ- 
ence probability density APdq S corresponding to the dominant pre-edge peak at P=6.5 eV, 
we can see that the largest effects purely connected with x-ray absorption actually occur at 
the nearest sulphur atoms. It is thus clear that nearest sulphurs definitely have their role 
in generating the s\\xy pre-peak. In fact, we found that omitting six nearest S atoms from 
a large 135-atoms cluster changes the calculated XANES drastically (making the very con- 
cept of pre-edge region inapplicable). Only the secondary-in-importance pre-edge peak at 
E ~7. 5-8.0 eV seems to be generated by scattering off second-shell Ti atoms, as it follows 
from AP^s curves in Fig. |l|. It seems, therefore, that the physical picture of the process 
which is responsible for the e\\xy pre-peak should consider the nearest sulphur atoms into 
account. Namely, although the total amount of the photoelectron density is small 
on sulphurs as compared with titaniums, its relatively modest variations with energy and 
polarization have a big impact on the XANES. 

Let us emphasize that, strictly speaking, this kind of analysis only informs about sen- 
sitivity of the photoelectron density near individual atoms to the changes of polarization 
vector direction. It does not testify about the physical mechanism which may stand behind 
the creation of the pre-peak electron states. So our conclusions do not in fact contradict the 
suggestions that those states arise due to hybridization of central Ti 4p and next-neighboring 
Ti 3d orbitals§ — both views may be rather complementary than opposing. 

Finally, let us note that the new look on the role of nearest sulphurs we take here may 
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be relevant to other non-centrosymmetric systems with a distinct XANES pre-peak, such as 
a colossal magnetoresistence material Lai-^Ca^MnOa (Ref. ^). 

B. Ti K edge XANES of rutile Ti0 2 

The Ti K edge XANES of rutile Ti0 2 was studied very intensively due to its interesting 



pre-edge structure and the debate still does not seem to be settled. BEiTE3 In this paper, 
however, we want to concentrate on the extended XANES up to ~120 eV, as a comprehensive 
study of the effect of individual atoms on spectral peaks in this region was performed.iH! 
In rutile TiC>2, the Ti atoms are located at the center of a distorted octahedron. The plane 
containing four equatorial oxygens O eq is parallel to the crystallographic c axis while two 
axial oxygens O ax lie inside the plane defined by axes a and b. The third coordination 
sphere is formed by two titanium atoms Ti± c , shifted from the central Ti by ±c. Instructive 
depictions of rutile Ti0 2 structure can be found, e.g., in figure 3 of Ref. 132], figure 2 of Ref. 



|35| , or figure 3 of Ref. |36| . 

Jeanne- Rose and Poumelleci! studied how the polarized Ti K edge XANES of TiC>2, 
calculated for a cluster of 25 atoms, is affected by variations in the positions of O ax and 
O eq atoms. That makes it possible for us to explore here to what extent does sensitivity 
of a spectral peak to a position of a certain atom imply high PEPD around that atom and 
vice versa. In accordance with the discussion in Sec. [111 A| , we will concentrate more on the 
difference probability APdos than on P^) itself. 



14 




FIG. 2. Theoretical polarized Ti K edge XANES of TiO"2 and corresponding PEPD curves. 
The meaning of the curves and symbols is analogous as in Fig. |l| 

In Fig. 0, theoretical polarized Ti K edge XANES of T1O2 is shown together with cor- 
responding pti\ AP^l, and APdos curves for the central Ti and its nearest neighbors. In 
order to connect with Jeanne- Rose and Poumellec,!! we consider a cluster of 25 atoms. The 
notation of spectral peaks C, D, E\, E 2 and E 3 is taken from Ref. 35, too. The ad-hoc 
Fermi energy would be around 10 eV in the energy scale of Fig. |2|, meaning that the wild 
oscillations at the beginning of the PEPD curves actually fall predominantly into the region 
of occupied states. 

Let us now compare the conclusions of Jeanne-Rose and Poumelleci! with the picture 
offered by the PEPD analysis. The equatorial oxygen O eq was found to influence both the C 
and the D peaks for both polarizations, the effect being significantly stronger on D than on 
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C and for the e\\c polarization than for the e\\a polarization.^ A counterpart to this effect 
can be found in the PEPD: A distinct peak in APdq S around O eq can be found at D and a 
mild feature at C for both polarizations; the D peak in AP^os is higher for the s\\c than for 
the e\\a polarization. 

When the position of the axial oxygen O ax varies, the calculated e\\a XANES is drastically 
altered at the D peak and not so much at the C peak, while the e\\c spectrum is changed 
at the C peak only.il A brief look at the O ax curves in Fig. ^ reveals that the largest 
APqos around the O ax atom is at the D peak for the s\\a spectrum and at the C peak for 
the e||c case. So indeed, the energy at which a particular atom affects the photoelectron 
probability density most prominently coincides in this case with the energy at which the 
XANES is significantly changed when that atom is moved. The analogy is nevertheless only 
a qualitative one — larger AP^os does not necessarily imply bigger changes of XANES peaks 
when O ax is moved. For example, the C resonance in APdq S for e\\c is stronger than the D 
resonance for e\\a, and yet the exact position of the O ax atom affects more the D spectral 
peak at e||c than the C peak at e\\a polarization.il On the other hand, the total probability 
densityP^ is larger at the D peak energy than at the C peak for both oxygens, reflecting 
correctly the greater sensitivity of the D peak to their positions. 

So the following intuitively plausible picture emerges: For comparing XANES peaks at 
different energies, the total probability pv> may be a good indication of their respective 
sensitivity to atomic positions. For investigating polarization effects, one has to resort to 
A P (j) 

i-i-l DOS' 

Unlike for the C and D peaks, a correspondence between XANES and APdos peaks 
cannot be established for the Ej maxima. Jeanne-Rose and PoumellecH found that the E2 
peak in the e||c spectrum is affected quite a lot by the O cq movement while peaks E x and 
E 3 are left basically intact for either polarization - a property that does not seem to have a 
counterpart in Fig. 0. The sensibility of a XANES peak to a movement of a particular atom, 
hence, does not necessarily imply a high localization of the photoelectron around that atom 
— that can be found only by a proper PEPD calculation. 

It is worth noting that if single-scattering dominates in generating a particular XANES 
peak (such as C for the e\\a polarization) ,E3 one can observe a distinct resonance in APq'qs &t 
corresponding energy. On the other hand, such a correspondence appears to be blurred for 
spectral peaks where a significant contribution from multiple-scattering is suspected (such 
as C peak at the Oa^atom for the e\\a polarization or the E2 peak at O eq and Ti± c for 
the e\\c polarization) x3 This seems to be plausible — a multiple-scattering nature of a peak 
emphasizes that it is generated not by a mere "presence" of the photoelectron near certain 
atom but rather by a complicated interference process. It remains to be explored to what 
extent this trend is a general one. 



IV. DISCUSSION 

Photoelectron probability density offers us a direct look on the spatial localization of 
those electron states which are seen by XAS. Contrary to a bit vague concepts like impor- 
tance of various atoms for the formation of a particular XANES peak or order of multiple- 
scattering which has to be accounted for, PEPD is a rigorously defined quantity with a 



16 



transparent physical interpretation. Hence even if there were no immediate practical appli- 
cations, PEPD analysis would still remain a valuable tool for understanding the physical 
processes which give rise to the XANES spectrum. Apart from this principal asset, there 
are two directions where PEPD analysis could contribute to solving concrete problems: in- 
vestigations of real structure of solids and investigations of their electronic structure. 

One of the obstacles to overcome when fitting the experimental XANES spectrum with 
calculated spectra of trial structures is that ones has, in general, a lot of atoms to move 
and hence it is difficult to identify those whose positions are most critical for the XANES 
shape. By providing a deeper insight into the "photoelectron trajectory", PEPD analysis 
may drop a hint for the most critical spots in advance. Intuitive arguments and practical 
experience (Sec. |T|) show that both and AP^os reflect the sensitivity of XANES to 
movement of individual atoms: Total probability density P^ is more relevant for comparing 
roles of different chemical species while APj&s 

is more indicative of the polarization and site 

dependence of XANES spectra. 

In electronic structure studies, the appealing feature of PEPD is that it directly inves- 
tigates those unoccupied states which are probed by XAS. When investigating the spatial 
localization of these states, one is thus not left relying on heuristic arguments like "which 
atom affects the spectrum most if moved" . One can even explore the angular-momentum 
character of XANES states by projecting the wave functions ^ } along Eqs. @-(g§. In 
that way, one can see not only from which atom a particular XANES peak arises but also 
from what type of orbital is comes from. 

The concept of PEPD relates, of course, not just to the XANES region but to the EXAFS 
part of the absorption spectrum as well. The reason why we mention only XANES here 
explicitly is that EXAFS oscillations can be properly analyzed with other tools. Note also 
that the PEPD analysis could be applied to photoelectron diffraction as well — one would 
just have to omitt the angular integration Jd 2 k in Eqs. (|7|), (|22|) and in related expressions. 



V. CONCLUSIONS 



Our analysis demonstrates that it makes sense to explore the probability that a pho- 
toelectron participating in a XANES process can be found at a specific place. The rele- 
vant quantity is the photoelectron probability density and it can be calculated as a sum 
of squares of wave functions describing elementary PED processes, weighted by normalized 
PED cross sections. When investigated as a function of the photoelectron energy, it exhibits 
a resonance-like structure and depends on the atom around which it is evaluated. The bulk 
of PEPD is dominanted by DOS-like effects, meaning that, e.g., hardly any polarization de- 
pendence can be seen in it unless the DOS-related portion is subtracted. In many cases, high 
difference probability AP^os around an atom may serve as an indication of high sensitivity 
of XANES towards the position of that atom for a given energy. The fine structure in PEPD 
does not copy the corresponding XANES: For some atoms and/or peaks, a visual connec- 
tion between XANES and PEPD can be established, for other features such a discernible 
connection is absent. This may be a manifestation of the interference nature of XAFS 
- the correspondence between XAS and PEPD peaks is more often observed for features 
which arise from single scattering than for multiple-scattering peaks. The spatial depen- 
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dence of PEPD thus provides information which is not equivalent to what can be learned 
from comparing theoretical spectra for various trial structures but rather is complementary 
to it. 

By performing a PEPD analysis for a Ti K edge of TiS2, we found that, contrary to earlier 
interpretations, the sulphur atom nearest to the absorbing titanium participates significantly 
in formation of the distinct polarization-dependence of the pre-peak. 
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APPENDIX A: OTHER EXPRESSIONS 
FOR PED AND XAS CROSS SECTIONS 



In this appendix, some other formulations of the equations presented in Sec. II B 



are 



shown, in order to facilitate connection with notations and definitions in other works. 
The photoelectron diffraction cross section, evaluated in Eq. ( f[8|) , can be expressed as 



da 



1 



dfifc 47T 



auk 



LLr 



i 2 2 



LLr 



(Al) 



By introducing the scattering matrix W of Eq. ( |14|) directly into Eq. ([18|) , one obtains 

2 



da 



E E W°AJ''Y* lf (k)e^ Rj -^D 

L jL" 



LL, 



+ 
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2 2 

-a uj 



EE<i"i £ "^»(fc)e ifc(RJ - R0) Ql 

L jL" 



LLr 



(A2) 



which is analogous to equation (22a) of Natoli et al£3 

As a lot of useful relations about XAS and PED within the wave function formalism can 



be found in Ref. [L6| we quote here the relation between the "outgoing" quantities used in 
that paper and the "incoming" quantities employed here. The scattering amplitude (3$ (L") 
of Eq. flT3"| ) is connected with the amplitude B L ^(L") of equation (2.28) of Natoli et aJlB via 

^\L") = (-l) m [B^",-m")Y(-l) m " (A3) 
and the scattering matrix W relates to the matrix 
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yij — 



I ■ H 



LL> 



of equation (3.3) of Natoli et atil as 

f-ir 



and to the r matrix of Refs. |37|,|38 



as 



'-D 



7«J 

^ {l-m){l' ,-m') 



ji 

T (l' ,-m')(£,-m) 



(A4) 



(A5) 



Employing the Z and r matrices, the x-ray absorption cross section of Eq. (^) can be 
written as 



a XAS = Airauk Im (pi L Zf v D VLc 

LL> 

+ -^a Vim (QI l ZTl>Ql>l c ) } 



(A6) 



or as 



-4naukY,{l™{DLLA L'D* m 

LL> 

+ -^a Vim (Q ll AI'Q*l>l c ) } 



(A7) 



Eqs. ( |A6|) -( [A7|) can be arrived at either by invoking the Green function formalism from the 
beginning or by applying the optical theorem on Eq. (|21|), as outlined by Natoli et a/.0 



APPENDIX B: 
PEPD IN CASE OF A SINGLE ATOM 



In calculating the single-atom probability density P a to(^), one can proceed just as in 
The only difficulty arises in expanding the solution of a single-atom Lippman- 



Sec. II B 



Schwinger equation inside (now) empty spheres around sites R ] . One cannot mechanically 
employ equations (p^)-(|I~4T) to find the amplitudes f3 ( f\k) } as the phaseshifts at non-central 
atomic sites are zero now. Instead, by applying formula for re-expanding ji(k\R p \) Yl(RP) 
and ni{k\R p \) Yl(RP) around different origins, one can arrive at the following expression for 



j e ,(k\R p \) cos 5^,, 



-n t ,{k\R p \)sm5$] Y V (R P )C^ L , . 



(o) 



(Bl) 



which ought to be inserted into Eq. (|22|) with the help of (0), together with replacement of 
the single-center radial wave function TZ^\kr) with its free-electron counterpart je(kr). 
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APPENDIX C: TECHNICAL DETAILS 
OF PEPD CALCULATIONS 



When evaluating the PEPD defined by Eqs. (p2|) and (25), we rely on the expression 



P ti) 



16r« ^fc^L E | E f'^) ( r)^„ ( fc) 



XAS 



H 



x / drr* 
lo 



k L L" 

Rf{kr) 



(CI) 



where the photoelectron diffraction cross section is evaluated from Eq. ( |T8| ) and the am- 
plitudes f3^'(L") from Eq. (|T3|). The angular fc-integration in ( |U1| ) could be performed 
analytically in principle. However, that would lead to such a proliferation of slowly conver- 
gent sums over angular momenta, that it is actually computationally more convenient to 
keep the fc-integral in the expression (|C1|) and to evaluate it numerically. 

The angular momentum sum J2l m Eq. (|C1|) descends from Eq. (|T0|) and converges quite 
quickly (it corresponds to a multi-center expansion in the terminology of Durham et at) £3 In 
all the cases analyzed in Sec. [Til], it was sufficient to cut this sum at £ max = 3. On the other 
hand, the sum which descends into Eq. ( |Cl"l ) from the expansion (pT3|), corresponds to a 
single-center expansion and converges only slowly. In this work it was cut at i- inc = 30, which 
we found to be a safe value. Taking £- mc = 20 would still lead to an acceptable accuracy, 
while increasing £ inc up to 40 did not induce significant changes with respect to the l- inc = 30 
case. 

The angular integration was performed in two "perpendicular" spherical coordinates 9, 

d> as 



d 2 k 



d9 sinO 

o Jo 



2n 



d<j) 



(G2) 



In Sec. |TJ, fifty points were used for numerical integration over the 9 coordinate and 
2 x 50 x sin 9 points for integration over 0. We checked that such a grid is dense enough to 
guarantee a sufficient accuracy. 
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